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Abstract. A metamorphic robotic system is an aggregate of homogeneous robot units which can 
individually and selectively locomote in such a way as to change the global shape of the system. We 
introduce a mathematical framework for defining and analyzing general metamorphic robots. This formal 
structure, combined with ideas from geometric group theory, leads to a natural extension of a configuration 
space for metamorphic robots — the state complex — which is especially adapted to parallelization. We 
present an algorithm for optimizing reconfiguration sequences with respect to elapsed time. A universal 
geometric property of state complexes — non-positive curvature — is the key to proving convergence to 
the globally time-optimal solution. 



1. Introduction. In recent years, several groups in the robotics community have been 
modeling and building reconfigurable or, more specifically, metamorphic robots (e.g., 
El El 122 ED- Such a system consists of multiple identical robotic cells in an underlying 
lattice structure which can disconnect/reconnect with adjacent neighbors, and slide, pivot, 
or otherwise locomote to neighboring lattice points following prescribed rules: see Fig. 11.11 
There are as many models for such robots as there are researchers in the sub-field: 2-d 
and 3-d lattices; hexagonal, square, and dodecahedral cells; pivoting or sliding motion: see, 
e.g., El IHHSl El El El 122 ESI and the references therein. The common feature of 
these robots is an aggregate of lattice-based cells having prescribed local transitions from 
one shape to another. 
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Fig. 1.1. Metamorphic systems may be built on a variety of lattice structures with sliding or pivoting 
motion. 



The primary challenge for such systems is shape-planning: how to move from one shape 
to another via legal moves. One centralized approach [HI El is to build a transition graph 
whose vertices are the various shapes and whose edges are elementary legal moves from one 
shape to the next. It is easily demonstrated that the size of this graph is exponential in the 
number of cells. 

We propose to extend the notion of a configuration space to metamorphic robots in a novel 
manner. The idea: consider the transition graph described above as a one-dimensional 
skeleton of a higher-dimensional cubical complex, the state complex. Assume that from a 
given state there are two legal moves which are physically independent (or, more suggestively, 
"commutative"): i.e., these moves can be executed simultaneously. In the transition graph, 
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this corresponds to the four edges of a square. For any pair of commutative moves, fill in 
the four edges of the graph with an abstract square 2-cell. Continue inductively adding k- 
dimensional cubes corresponding to fc-tuples of physically independent motions. The result 
is a cubical complex which has several advantages over the transition graph: 

1. Simplicity. The state complex is often simpler than the transition graph: e.g., 
the 1-d graph of an n-dimensional cube has n2 n ~ 1 edges. This figure belies the 
simplicity of the single cube. 

2. Speed. Geodesies on this complex cut across the diagonals of cubes whenever 
possible. One performs all possible commutative motions simultaneously, maxi- 
mizing parallclization and yielding a speed-up by a factor equal to the number of 
coordinated motions. 

3. Shape. The global geometry /topology of the state complex carries information 
about the metamorphic system. For certain examples, the topology of the state 
complex "converges" upon refining the lattice. In addition, only special geometries 
can be realized as the state complex of a local metamorphic system: commutativity 
in reconfiguration leads to an abhorrence of positive curvature in the state complex. 

Sections |21 through 01 give definitions and examples of [abstract] metamorphic systems and 
their state complexes. The next two sections (Sections EE} detail topological and geometric 
features of the state complex. 

For large systems, the problem of computing the state complex and designing geodesies in 
order to perform shape planning is computationally infcasible, primarily because the size 
of the complex is often exponential in the number of robot cells. In addition, any control 
scheme induced by geodesic construction is necessarily centralized. Several researchers have 
begun building decentralized control algorithms for shape planning jSJ 1211 121)1 1241 125| . Such 
algorithms have the advantage of speed and scalability; however, the reconfiguration paths 
are typically not optimal. 

As an application of our techniques, we present in Sections 17181 an algorithm for trajectory 
optimization which takes as its argument an arbitrary edge path in the transition graph. 
Algorithm 18 . 21 then performs a type of curve shortening within the state complex. A deep 
theorem about the curvature of all state complexes (Theorem l6.2|l is then used to prove that 
this algorithm returns a shape trajectory which is the global minimum obtainable from this 
path with respect to elapsed time. 

Our definitions and theorems are phrased for systems involving "discrete" reconfiguration. 
More general types of robots which employ continuous reconfiguration for locomotive gaits 
(such as the Polybot developed by M. Yim's lab at Xerox PARC) are not covered by our 
definitions. We note, however, that certain locomotive reconfigurable robots can be thought 
of as lattice-based tiles by amalgamating subsystems |T3]. In addition, our definitions can 
easily be extended to more general non-lattice reconfigurable systems jHJ. 

2. A mathematical definition. While it is easy to generate examples of what is 
meant by a metamorphic robot, it is more challenging to write a clean mathematical defi- 
nition. We propose a set of definitions which is broad enough to include some non-obvious 
examples. The paper [H] suggests a similar type of structure using cellular automata rule 
sets. 

A local metamorphic system is a collection of states on a lattice, where each state is thought 
of as an indicator function for the aggregate. Any state can be modified by local rearrange- 
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ments, these local changes being coordinated by a catalogue of models realized under the 
actions of isometries into the workspace. The adjective "local" refers to legality criteria: 
anywhere in the workspace at which a local change from the catalogue can be applied, it 
is legal to do so. To incorporate obstacles and basepoints into our systems, we distinguish 
between the amount of information needed to determine the legality of an elementary move 
(the "support" of the move) and the precise place in which modules are actually in motion 
(the "trace" of the move). 

Definition 2.1. Let C denote a lattice in M. k and let W C C be some workspace. The 
catalogue C for a local metamorphic system on W is a collection of generators. Each 
generator <fi £ C consists of (1) the support, SUP(0) C C; (2) the trace of the move, 
tr(0) C SUp(0); and (3) an unordered pair of local states Uqj : SUP(0) — > {0,1} satisfying 1 



(2.1) U 



= Ui 

SUP(0)-TR(0) 



SUP(0)-TR(<£) 

Otherwise said, the local states are equal on SUP(0) — TR(</>). 

Definition 2.2. An action of a generator <\> G C is a rigid translation €> : SUP(0) <— > W. 
Given a state U : W — > {0, 1}, an action $ is said to be admissible at U if Uq — U o $ . 
In this case, we write 



U : on W - *(SUP(0)) 
Ui o <I> 1 : on <&(sup(0)) 



(Note that & is left out of the notation.) 

Definition 2.3. A local metamorphic system on W is a collection of finite states 
{U a : W — > {0, 1}} closed under all possible admissible actions of generators in the catalogue 
C. (A state is finite if U = 1 at only finitely many points of W.) 

To repeat, the catalogue and the workspace are the "seeds" for a local metamorphic system. 
From this pair, all possible translations of the supports into W yield the actions. Then, a 
collection of states on the workspace is a local metamorphic system if, whenever an action 
<fr of a generator ^ on a state U is admissible, then the corresponding state <fr[U] is also 
included. 

A metamorphic system with obstacles O C W satisfies in addition 
(2.2) *(tr($) n O = 0, 

for each cf> e C. Obstacle sets count as legal positions for determining the admissibility of 
a move (the support of an action may intersect 0), but no motion of metamorphic agents 
may incorporate the obstacle sites (the trace of an action must not intersect O). 



3. Examples. Some of the following examples are inspired by metamorphic robots 
already developed; other examples are more abstract. 

Example 3.1. [2-d hex with pivots] We present two slightly different catalogues, each with 
six generators (or one, up to discrete rotations), in Fig. 13. II Both of these systems, modeled 
after that of jxj, have local moves which pivot a planar hexagon about a neighbor. For all 
generators presented, the trace is equal to the two central hexagons. In the first system, 



All generators are assumed to be nondegenerate in the sense that Uo ^ Ui. 
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the support is chosen so that the aggregate does not change its topology but only its shape. 
The slightly smaller support of the second catalogue allows for local topology changes. To 
model a fixed "base" cell (which is, say, affixed to a power source as in [TJ\), one establishes 
this cell as an obstacle O. 

Fig. 3.1. Two different catalogues for a 2-d hexagonal lattice system with pivots. Black cells are 
occupied, white are unoccupied. The local states Uq and Ui are shown for each generator. 

Example 3.2. [2-d square lattice] In Fig. 13.21 we display a generator for a planar system 
in which rows [as pictured] and columns [not pictured] of an aggregate of square cells can 
slide. There are in fact several generators represented in "shorthand," one for each k > 0. 
A dot inside a cell indicates that it can be either occupied or unoccupied, but if occupied, 
then its neighbor (indicated by an arrow) must also be turned on. This condition guarantees 
that the aggregate does not disconnect (even locally) under slides. The trace of this set of 
generators is the entire middle row except the two endpoints. To keep the catalogue finite, 
one would include only those generators with k < N, where N is the number of occupied 
cells in any state. 
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Fig. 3.2. The row-sliding generators for a sliding- squares system. 



Example 3.3. [2-d articulated planar arm] Consider as a workspace W the set of edges 
in the planar integer lattice. The catalogue consists of two generators, pictured in Fig. 13.31 
Beginning with a state having N vertical edges end-to-end, the metamorphic system thus 
generated models the position of an articulated robotic arm with fixed base which can (1) 
rotate at the top end and (2) flip corners as per the diagram. 



Fig. 3.3. A positive articulated robot arm [left]. One generator [center] flips corners and has as its 
trace the central four edges. The other generator [right] rotates the end of the arm, and has trace equal to 
the two activated edges. 

If one includes rotations of these generators, more intricate types of configurations are possi- 



COMPLEXES FOR METAMORPHIC ROBOTS 



5 



ble, including deadlocked configurations. Highly self-reconfigurable examples with multiple 
interacting arms can be realized by adding new generators: an "attach-detach" generator 
which allows endpoints of arms to merge (thus yielding a "marked point" at the attach- 
ment); and a "sliding" generator which allows this marked point to slide, having the effect 
of allowing the attached arms to trade segments. These generators are pictured (up to 
Euclidean symmetries) in Figure l3~4l 




Fig. 3.4. Additional generators allow for attachment and detachment of arm endpoints [left, center] 
and sliding of the coupling points [right]. 

One of the benefits of writing down a rigorous definition of a metamorphic system is the 
discovery of systems which have little resemblance to the systems of, say, Fig. 11.11 In 
particular, our definitions easily extend to metamorphic systems which are not lattice-based; 
the following example is especially interesting. 

Example 3.4. Consider a finite graph T in which every edge is assigned a length of one. 
(Every graph can be embedded in some K fc so as to have this property.) The catalogue 
consists of a single generator whose support and trace are precisely the closure of a single 
"abstract edge." The local states of this generator consist of the pair Uq and Ui which 
evaluate to 1 on one of the endpoints and on the other. The actions in this case are 
length-preserving maps from the abstract edge into T. The metamorphic system generated 
from a state U on T with N vertices evaluating to 1 mimics an ensemble of N unlabeled 
non-colliding Automated Guided Vehicles on T, cf. 

More abstract examples of metamorphic examples include spaces of triangulations of poly- 
gons with edge-flipping as the generator, examples arising from word representations in 
group theory, and certain multi-step assembly processes [3] . 

4. The state complex. In the robotics literature, one often models a configuration 
space for a metamorphic system with a transition graph which represents actions of elemen- 
tary moves on states. That is, the vertex set is the collection of all states {U a }, and the 
edges are unoriented pairs of states which differ by the action of one generator. Transition 
graphs are discussed for shape-planning in several particular cases in the literature (planar 
hex case: jSJEl)- Our departure is to make the transition graph the 1-skeleton of a cubi- 
cal complex (an analogue of a simplicial complex, but made out of abstract cubes) which 
coordinates parallel or "commutative" motions. 

Definition 4.1. In a local metamorphic system, a collection of actions of (not necessarily 
distinct) generators {(^a*)^a*)} * s sa id t° commute if 



(4.1) 



*a s (tr(0 Qi )) n * Qj (sup(<?(> ttj )) = V^ + j. 
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Example 4.2. Two simple examples suffice to illustrate the difference between commuting 
and noncommuting actions. First, consider the pair of commuting moves for a planar 
hexagonal pivoting system, as represented in Fig. 14.11 [left]. Compare this with a planar 

i i I T 




Fig. 4.1. Examples of commuting [left] and noncommuting [right] actions in planar systems. 

sliding block example as illustrated in Fig. 14. II [right]. Although the pair of moves illustrated 
forms a square 2 in the transition graph, this particular pair of actions does not commute. 
Physically, it is obvious why these moves are not independent: sliding the column part- way 
obstructs sliding a transverse row. Mathematically, this is captured by the traces of the 
actions intersecting. 

The state complex has an abstract fc-cube for each collection of k admissible commuting 
actions: 

Definition 4.3. The state complex S of a local metamorphic system is the following 
abstract cubical complex. Each abstract k-cube e^ k ' of S is an equivalence class [U ; {^ ai )^ = i\ 
where 

1. is a k-tuple of commuting actions of generators cf> ai ; 

2. U is some state for which all the actions (^& ai )k =1 are admissible; and 

3. [Uo; (&ai)i=i\ = [Ui', (*/3*)i=i] */ an d on ly if the list (ft) is a permutation of {ai) 
and U — U\ on the set W — \J i <fr ai (sup(0 ai )) . 

The boundary of each abstract k-cube is the collection of 2k faces obtained by deleting the 
i th action from the list and using U and <j) ai [U] as the ambient states. Specifically, 

k 

(4.2) 0[CM*ai)iU] = U ([^(^WduI^jC/];^)^]). 

i=l 

It follows easily that the fc-cells are well-defined with respect to admissibility of actions. The 
proof of the following obvious lemma is given in detail to flesh out the previous definition. 

Lemma 4.4. (a) The O-dimensional skeleton of S, S^ ', is the set of states in the reconfig- 
urable system, (b) The 1-dimensional skeleton of S, is precisely the transition graph. 



2 The individual robotic cells are not labeled: only the shape of the aggregate is recorded. 



COMPLEXES FOR METAMORPHIC ROBOTS 



7 



PROOF: (a) Vertices of S consist of equivalence classes consisting of zero (i.e., no) actions 
of generators up to permutation, together with a state defined on the complement of the 
supports of the actions. As there are no actions, each 0-cell is precisely a single state of the 
reconfigurable system. 

(b) A 1-cell of S is an equivalence class of the form [U ; (<&)] . The only other representative of 
the equivalence class is (*)]; hence, the 1-cells are precisely the edges in the transition 

graph. Clearly, the boundary of [U; (<&)] is the pair of 0-cells [U; (■)] and [0 [J7] ; (•)]. o 

For small numbers of cells, it is easy to illustrate the state complex. 




Q 



V 

Fig. 4.2. The planar sliding square example [left] with two movable blocks [in black] and a p-by-q 
obstacle set [in grey] yields a state complex S that is topologically a circle [right]. 

Example 4.5. Consider the 2-d square lattice row/column sliding system whose catalogue 
is illustrated in Fig. 13.21 If we consider a system with (occupied) obstacles in the form of a 
p-by-q rectangle generated from the state of Fig. 14.21 [left], one obtains a planar transition 
graph with A(pq + 1) + 2(p + q) vertices and 8(pq + 1) — 2(p + q) edges. In contrast, the 
state complex is that of Fig. 14.21 [right]: this is topologically a circle, corresponding to the 
fact that the pair of free squares can circulate about the obstacle set through a sequence 
of slides. The large 2-d regions correspond to states in which the two free squares are on 
separate (but adjacent) sides of the obstacle set. 

Example 4.6. Consider the planar hex system of Fig. 13.11 [left] with a workspace W 
consisting of a long channel of four rows, the bottom row being filled obstacles, and the top 
row being non-occupied obstacles. A line of cells in one of the two free rows can "climb 
across" one by one, yielding the state complex of Fig. I4.3| right] (cf. the algorithm of |2(J|). 
This state complex is contractible for any length column. 




Fig. 4.3. For a line of hexagons filing out of a constrained tunnel [left], the state complex is contractible 
[right]. Black cells represent occupied obstacle sets; the top row is an unoccupied obstacle set. A fairly short 
tunnel is shown: for longer tunnels, the state complex is still contractible, but of higher dimension. 
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Example 4.7. The state complex associated to the positive articulated robot arm of 
Example 13.31 in the case N — 5 is given in Fig. 14.41 Note that there can be at most 
three independent motions (when the arm is in a "staircase" configuration); hence the state 
complex has top dimension three. Notice also that although the transition graph for this 
system is complicated, the state complex itself is topologically trivial (contractible). 




Fig. 4.4. The state complex of a 5-link positive arm has one cell of dimension three, along with several 
cells of lower dimension. 

Example 4.8. In the system of Example 13.41 with the graph being a K§ (the complete 
graph on five vertices) and N = 2, the state complex is a two-dimensional closed surface. 
A simple combinatorial argument (as in QUI]) reveals that the Euler characteristic is —5, 
implying that the state complex is non-orientable. If the AGV's are labeled, the state 
complex becomes a closed orientable surface of genus 6. 

5. The topology of S. If one looks at a transition graph without knowing the particu- 
lars of the metamorphic system, very little information can be extracted. This paper argues 
that completing the transition graph to the state complex is "natural" — the state complex 
simplifies the transition graph and endows it with topological and geometric content. 

Our first example of naturality is motivated by the desire to build metamorphic systems with 
large numbers of micro- or nano-scale cells. While large numbers of cells would yield a type 
of continuum-limit convergence on the dynamics of shape change, the resulting transition 
graphs have no such convergence. The size of the transition graph goes up exponentially in 
the number of cells; more ominous is that the topology of the transition graph (the number 
of basic cycles) blows up as well. This is not always so with the state complex: in certain 
key examples, the topology of S is either invariant or converges to a limiting type. 

We have already seen one such example of this stabilization. In Example 14.51 the state 
complex of a pair of squares sliding along a rectangular obstacle of size p-hy-q is topologically 
a circle, independent of p and q. This can be interpreted as a type of convergence: consider 
the effect of refining the underlying lattice structure, increasing p and q while maintaining a 
pair of sliding squares. Then the dimension of the state complex remains the same, as does 
the topological type of the space. Intuitively speaking, the "limit" as this refining process 
is repeated yields a "topological" configuration space of two points sliding smoothly along 
the boundary of a rectangle, which can dock or un-dock at the corners. 

Example 5.1. Recall the state complex associated to the metamorphic system of N points 
on a graph T, Examplc l3.4l Consider a refinement of T which inserts additional vertices along 
edges. It follows from the techniques of pQ that the state complex of this refined system 
has the same topological type (up to homotopy equivalence) after a fixed bound on the 
refinement (N additional vertices per edge). Furthermore, this "stabilized" state complex 
is in fact homotopic to the topological configuration space of N non-colliding points of T — 



COMPLEXES FOR METAMORPHIC ROBOTS 



9 



precisely what one expects as the number of refinements goes to infinity. 

Example 5.2. Recall the positive articulated robot arm of Example 13. 31 Consider a refine- 
ment of the underlying lattice which shrinks the lattice by a factor of two (or, equivalently, 
which inserts an additional joint in the middle of each edge). This is a more dramatic 
change since the dimension of the state complex doubles. Nevertheless, the topological type 
is invariant: the state complex remains contractible. 

Proposition 5.3. Let Sn denote the state complex of the positive articulated arm from 
Examvle W.'ft with N segments. The complexes Sn are all contractible. 

PROOF: For these articulated arms, there is a nice inductive structure on the state 
complexes. Fixing N, each state (vertex) in Sn is represented as a length N word in the 
symbols x and y, where x denotes the arm going to the right and y denotes the arm going 
up. In this language, the two generators are (1) transposing a subword xy <-> yx, and (2) 
changing the last letter of the word. 

Consider the subcomplex X C Sn consisting of all cells whose vertices have words beginning 
with the letter x. Likewise, let Y denote the subcomplex all of whose vertices begin with the 
letter y. These subcomplexes are each a copy of Sn-i which we may assume inductively is 
contractible. One passes between the subcomplexes X and Y only when a move exchanges 
the initial two letters of the word from xy to yx. The connecting set is thus homeomorphic to 
Sn-2 x [—1,1] and attached to X and Y along Sn-2 x { — 1} an d Sn-2 x {1} respectively. 
Again, by induction, these sets are contractible. A pair of contractible sets joined along 
contractible subsets is contractible. o 

This should come as no surprise: in the limit as N — > oo, the reconfigurable system approxi- 
mates the configuration space of a smooth curve of fixed length which is positive in the sense 
that the curve is always nondecreasing in the horizontal and vertical components. That the 
(infinite dimensional) space of such smooth curves is contractible is easily demonstrated: 
given any such curve with endpoint fixed at the origin in the plane, pull the other end along 
the straight line connecting it to the origin until the strand is taut. Then, rotate the line 
segment rigidly until it is, say, vertical. This is a continuous deformation on the space of all 
smooth positive curves of fixed length to a single vertical segment. 

It is certainly not the case that an arbitrary reconfigurable system possesses such convergence 
properties: the manner in which one refines the states is important. Still, we conjecture 
that state complexes can often be viewed as "discretizations" of some underlying smooth 
configuration space. 

6. The geometry of S. There are several natural ways to measure distances in state 
complexes. We first discuss the geometry arising from considering each cube of S to be 
Euclidean (i.e., flat), with unit side length; we call iS with this metric a Euclidean cube 
complex. However, this does not imply that the complex, as a whole, is flat. Indeed, non- 
zero curvature can be concentrated at places where several cells meet. A simple example 
appears in Fig. 16.11 here, a surface built from flat 2-cells can be seen to have curvature 
which depends on the number of 2-cells incident to a vertex. Four incident cells implies 
zero curvature; three cells implies positive curvature; and five or more cells implies negative 
curvature. For a two-dimensional complex, this is equivalent to computing the total angle 
about a vertex. 

Such an extension of curvature to general metric spaces is made precise in Gromov's work 
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Fig. 6.1. Curvature about a vertex in a 2-d Euclidean cubical complex: positive, zero, and negative. 

on curved metric spaces |12| (extending the classical work of Alexandrov, Busemann, and 
others) in which triangles with geodesic edges are used to measure curvature bounds. In 
brief, let X be a metric space and p € X a point. To bound the curvature of X at p, consider 
a small triangle T about p with geodesic edges of length a, b, and c. Build a comparison 
triangle T' in the Euclidean plane whose sides also have length a, b, and c respectively. 
Choose a geodesic chord of T and measure its length d. In T", measure the length d! of the 
chord whose endpoints correspond to those of the chord in X . 

Definition 6.1. A metric space X is nonpositively curved (or NPC,) if for every 
sufficiently small geodesic triangle T and for every chord of T , it follows that d < d' . 




Fig. 6.2. Comparison triangles measure curvature bounds. 

In other words, geodesic chords are no longer than Euclidean comparison chords. It should 
be stressed that the NPC property is very special and highly desirable. Indeed, being NPC 
implies a variety of topological consequences reminiscent of smooth nonpositively curved 
manifolds. 

Despite the variety of (local) metamorphic systems, all state complexes share this special 
geometric property. 

Theorem 6.2. The state complex S of any local metamorphic system is nonpositively 
curved. 

The proof of this theorem is simple, but requires some additional machinery 

Definition 6.3. Let X denote a complex (either simplicial or cubical) and let v denote 
a vertex of X. The link of v, £k[v], is defined to be the abstract complex which has 
one fc-dimensional simplex for each (k + l)-dimensional cube in X incident to v. The 
boundary relations are those inherited from X: namely, the boundary of a fc-simplex in 
£k[v] represented by a (fc + l)-cube in X is the set of all simplices represented by the faces 
of the (k + l)-cube. 

Links can be thought of as a simplicial version of the locus of points a small fixed distance 
from the vertex v. 



Definition 6.4. A Euclidean cube complex X satisfies the link condition if, for each 
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vertex v G X, £k[v] satisfies the following: for each fc, if any k + 1 vertices in £k[v] are 
pairwise connected by edges in £fc[u], then those vertices bound a unique fc-simplex in £fc[t>]. 

An important and deep theorem of Gromov 12| asserts that a Euclidean cube complex is 
nonpositively curved if and only if it satisfies the link condition. This criterion makes it 
easy to prove Theorem 16. 21 

PROOF OF theorem 16.21 Let U denote a vertex of S. Consider the link £k[U]. The 
0-cells of the £k[U] correspond to all edges in S incident to U; that is, actions of generators 
admissible at the state U. A fc-cell of £k[U] is thus a commuting set of k + 1 of these actions 
based at U. The interpretation of the link condition for a state complex is as follows: if at 
U one has a set of k + 1 admissible actions, {(</> Qi , <&a)}i=i, of which each pair commutes, 
then the full set of k + 1 generators must commute [existence of the fc-simplex in the link] . 
Furthermore, they must commute in a unique manner [uniqueness of the fc-simplex in the 
link]. 

The proof of existence follows directly from Definition 14.11 any collection of pairwise com- 
mutative actions is totally commutative. We therefore have a (fc + l)-dimensional cell in 
S which is the equivalence class [U; (^aj^f 1 ]- This is the representative in S of the fc- 
simplex in £k[U]. To show uniqueness, consider any other fc-simplex in £k[U] which has the 
same vertex set. This must correspond to a (fc + l)-dimensional cube in S with actions 
{4> ai , &ai)i=i ( U P t° some permutation) based at U. From Defmition l4.3l this cell must be 
the same equivalence class, namely [U; (*a 4 )i=i]- * 

Example 6.5. To see how the NPC property can fail, consider the following non-local 
metamorphic system. Recall the generator for the planar hexagonal system presented in 
Fig-IOIright] which (along with its rotations) allows for local disconnection of the aggregate. 
If we add to this system a global rule that requires the aggregate to be connected (for, say, 
considerations of power transmission) , then we no longer have a local system, and positive 
curvature may exist. Fig. 16. 31 shows a configuration in which three actions of local generators 
act to disconnect the aggregate locally but not globally. These actions commute pairwise, 
and any two do not disconnect globally. However, performing all three actions disconnects 
the space and leads to an illegal state. Therefore, the corresponding state complex for this 
non-local system has a "corner" of positive curvature as a local factor. Note: the state 
complex will not be two-dimensional here, but will be locally a product of this with another 
cubical complex. The positive curvature persists. 



Fig. 6.3. A generator which allows for local disconnection [left] admits configurations [center] for 
which pairs of actions leave the aggregate connected, but triples do not; this non-local system has positive 
curvature in the state complex at this state [right]. 
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7. Geodesies and time on S. Besides displaying a unifying geometric feature, the 
nonpositive curvature has implications for path-planning, and hence, the shape-planning 
problem. 

Corollary 7.1. Each homotopy class of paths connecting two given points of a state 
complex contains a unique shortest path. 

PROOF: This is well-known for NPC spaces. The only difficulty lies in proving that the 
d < d! inequality for sufficiently small triangles implies the d < d! inequality for all geodesic 
triangles which are contractible in the space: see g|. Assume, then, that this inequality 
holds in general, and consider a pair of distinct homotopic shortest paths from points p to q 
in S. Choose any point r on one of the two paths. The two halves of this path bisected by r 
are themselves shortest paths from p to r and r to q respectively. Thus, there is a geodesic 
triangle in S whose comparison triangle in the Euclidean plane is a degenerate straight line 
from p' to q' . The d < d! inequality applied to the segment from r to its corresponding point 
on the other geodesic shows that these points coincide; therefore the two shortest paths from 
p to q in S are identical, o 

Fig. I7.1| left] gives a simple example of a 2-d cubical complex with positive curvature for 
which the above corollary fails. 

This corollary is a key ingredient in the applications of NPC geometry to path-planning on 
a configuration space, since one expects geodesies on S to coincide with optimal solutions 
to the shape planning problem. However, in the context of robotics applications, the goal of 
solving the shape-planning problem is not necessarily coincident with the geodesic problem 
on the state complex. Fig. I7.2[ lcft] illustrates the matter concisely. Consider a portion of a 
state complex S which is planar and two-dimensional. To get from point p to point q in S, 
any edge-path which is weakly monotone increasing in the horizontal and vertical directions 
is of minimal length in the transition graph. The true geodesic is, of course, the straight 
line, which is hardly a discrete object and thus more difficult to compute. 



Fig. 7.1. [left] A 2-d cube-complex (with positive curvature) possessing distinct homotopic shortest 
paths between a pair of marked points — note all edges have length one; [right] a 2-d cube-complex (with 
positive curvature) possessing an edge path (the thick line on the boundary) which is a locally (but not 
globally) shortest path. Any cube path near this path is strictly longer. 

Given the assumption that each elementary move can be executed at a uniform maximum 
rate, it is clear that the true geodesic on S is time-minimal in the sense that the elapsed 
time is minimal among all paths from p to q. However, there is an envelope of non-geodesic 
paths which are yet time- minimizing. Indeed, the true geodesic in Fie. l7.2| lcft] "slows down" 
some of the moves unnecessarily in order to maintain the constant slope. 

This leads us to dehne a second metric on 5, one which measures elapsed time. Namely, 
instead of the Euclidean metric on the cells of S, consider the space S with the norm on 
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Fig. 7.2. [left] The true geodesic lies within an envelope of time-minimizing paths. No minimal 
edge-paths are time- minimizing; [right] a normal cube-path in an NPC complex (after 1181 ). 

each cell. (This is also called the supremum norm: a vector is measured by the maximum 
of its components in each coordinate direction.) The geodesies in this geometry represent 
reconfiguration paths which are time minimizing. Using the results of |18| . one can prove 
that these geodesies are easily described using the notion of a cube path. 

Definition 7.2. A cube path from a vertex vq to a vertex vn is an ordered sequence of 
closed cubes {Ci}^ in S which satisfy (1) Ci Pi Ci+i = Vi for some vertex Vi; and (2) Ci 
is the smallest cell of S containing Vi-i and Vi. A cube path is said to be normal if in 
addition (3) for all i, St(Ci) D C;+i = Vi, where St(Cj) is the star of d (the union of all 
closed cubes, including Ci, which have Ci as a face). 

Roughly speaking, a normal cube path is one which uses the highest dimensional cubes as 
early as possible in the path. 

Theorem 7.3. In any state complex S, there exists a unique normal cube path from p to 
q in each homotopy class. The chain of diagonals through this cube path minimizes time 
among all paths in its homotopy class. 

This result follows from Theorem 16.21 and a result of JS| ( or > see the constructive algorithm 
of Section IS] below for a proof of existence). 

This result (or, alternatively, the constructive algorithm in the next section) also implies 
that, in the category of homotopic cube-paths, there is no such thing as a strictly local 
minimum for length: if a cube path in S has the property that all nearby cube paths are 
longer (more cubes), then this cube path is a normal cube path and indeed is the shortest 
cube path within the homotopy class of paths. See Fig. I7.1| right] for a simple example of a 
cube complex with some positive curvature that has a cube path which is locally of minimal 
length among cube paths, but not globally so. 

This is very significant. By beginning with any path in the state complex, say, one obtained 
via a fast distributed algorithm for shape-planning, one can employ a gradient-descent curve 
shortening on the level of cube paths. The above results imply that any algorithm which 
monotonically reduces length must converge to the shortest path (in that homotopy class) 
and cannot be hung up on a locally minimal cube path. The presence of local minima is a 
persistent problem in optimization schemes: nonpositive curvature is a handy antidote. 

8. Optimizing paths. In cases where the state complex is sparse — local principal 
cells being of high dimension with few neighbors — it is possible to compress the size 
of the transition graph significantly. However, shape-planning via constructing all of S 
and determining geodesies is, in general, computationally infeasible: the total size of the 
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state complex is often exponential in the number of cells in the aggregate. We therefore 
assume that some shape trajectory has been determined (by a perhaps ad hoc or distributed 
method) and turn to the problem of optimizing this trajectory. The nonpositive curvature 
of iS allows for a time-optimization which does not require explicit construction of S. We 
detail an algorithm for transforming any given edge-path to a time-optimal normal cube 
path. 

From Definition 14.31 an n-dimensional cube C of S can be represented by a set of n com- 
mutative actions along with an admissible state U. In the following algorithm, we 
suppress the state for notational convenience and consider C as the set of actions. Addition 
and subtraction is defined by adding or taking away admissible commutative actions to the 
list. 

Algorithm 8.1 (TimeGeodesic). Given: a cube path C = {Ci} i=1 ,, N in S. 

1: Let N := Length(C). 

2: Call ShrinkCubePath(C). 

3: If Length(C) < N then 1: else stop. 

ALGORITHM 8.2 (ShrinkCubePath). Given: a cube path C = {Ci}i=i..N in S. 

1: Let i = 1. 

2: Let X := Commute(Ci; Cj+i)- 

3: Update d := C t + X; 

4: Update C i+1 := C l+1 - X. 

5: Call ExciseTrivial(Ci+i). 

6: Call CommonEdge(Ci_i, C%). 

7: Call ExciseTrivial(Cj_i; Ci). 

8: If X = then % = i + 1. 

9: If Cj+i ^ then 2: else stop. 

Subroutine 8.3 (Commute). Given: a pair of cubes Cj = {$£ 4 }i=i..m ar| d Cj+i = 
l* i+1 l i 

1: Let5 : =Usup(*J.) 
2: LetT:=U i TR(*i 1 )_ ^ 
3: For i — h.n return <I?^ + if 
3.1: 5riTR(*^ j + 1 ) = ; and 

3.2: rnsup($|+ 1 ) = 0. 

Subroutine ExciseTrivial checks a cube for an empty list, removes these from the path, 
and reindexes the cube path, thus reducing the length. Subroutine CommonEdge checks a 
cube against adjacent cubes in the sequence for common edges and deletes these, returning 
a genuine cube path (recall Dcnnition l7.2fl . 

Subroutine Commute takes as its argument a pair of cubes Ci and Cj+i and returns those 
elements of C^+i which commute with all elements of Cj. The following lemma is a manip- 
ulation of the definitions. 

Lemma 8.4. The result o/Commute(Ci; Cj+i) is precisely the set of edges in St(Cj) PlCj+i. 

Theorem 8.5. Given a cube path C, Alaorithm \8.1\ computes a globally time-optimal path 
in the homotopy class ofC. 
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Fig. 8.1. Three rounds of shortening an edge path C. Cube path length stabilizes after two calls to 
ShrinkCubePath. Three more calls produces a normal cube path. 



PROOF: Recall that in this setting, the length of a cube path refers to the number of cubes 
in the path; this equals the time required to execute the reconfiguration. 

Repeated calls to Algorithm ShrinkCubePath eventually leave a cube path fixed. To show 
this, note that the integer-valued function f(C) := J^i • dim(Ci) decreases in each call to 
ShrinkCubePath which changes the path. Lemma f8.4l implies that ShrinkCubePath leaves 
a cube path fixed if and only if it is a normal cube path. From Theorem 17.31 this yields a 
time-optimal path. 

However, Algorithm TimeGeodesic calls ShrinkCubePath only until the length of the cube 
path is unchanged. Thus it remains to show that once the length of a cube path in unchanged 
by a call to ShrinkCubePath, this is equal to the length of the associated normal cube path. 
To show this, assume that ShrinkCubePath stops without excising any trivial cubes (i.e., 
shortening the length). We show that further calls merely redistribute actions earlier in 
the list C without deleting any permanently. Without a loss of generality, assume that the 
first call to ShrinkCubePath does not change the length of the cube path, but that in the 
subsequent call, the cube C^+i is eliminated via Ci commuting with a portion of Ci-\. see 
the schematic of Fig. 18.21 




Fig. 8.2. Schematic of exchanging actions in sequential rounds. Each dimension in this cartoon 
represents a multi- dimensional cube. 



Since C^+i is not eliminated in the first call, there must have been actions in Ci which 
did not commute with those of Ci-i but which are subsequently pushed backward in the 
second call. In addition, all the actions of C^+i are contained in ("parallel to" in the figure) 
the subset of Cj_x which commutes with Ci in the second call. This being the case, the 
commutativity of Ci+i with Ci would have occurred in the first call: contradiction, o 
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It is not surprising that Algorithm 18.11 converges to a locally time-minimal solution. The 
interesting implication is that the cube-path obtained is the global minimum for all paths 
obtainable from the initial: there is no way to get a quicker reconfiguration path by first 
lengthening the path and then shortening. This is the boon of nonpositive curvature. Fig. 18. II 
gives an example of paths generated by Algorithm 18. II 

The complexity of Algorithm l8.1l depends on both the catalogue of generators for the system 
and the length N of the input edge path. However, the dependence on the catalogue 
arises only in Subroutine Commute, which tests various pairs of subsets of the workspace for 
disjointness; the sizes of these subsets are governed by the generators. Assuming a fixed 
catalogue, then, there is a constant bound on the time required for a single call to Commute, 
and we may analyze the complexity of Algorithm 18 . II as a function of N alone. 

Notice that Algorithm 18 . II calls ShrinkCubePath at most N times, since it stops as soon as 
the length of the cube path fails to decrease. Each time an individual loop of ShrinkCubePath 
is executed, the quantity Y] dim(C 7 ) decreases. Since this quantity equals N initially, 
ShrinkCubePath calls the subroutines Commute, ExciseTrivial, and CommonEdge each at 
most N times. Each of these has a constant running time, with Commute being the only 
one depending on the particulars of the metamorphic system. Therefore the complexity of 
the entire Algorithm 18.11 is 0(N 2 ), with the constant determined by the running time of 
Commute. 

The worst case for Alaorithm lS. ll seems to be achieved in a totally flat 2-d state complex by a 
path consisting of two perpendicular segments of length N/2. (This case requires N 2 /2 runs 
of ShrinkCubePath.) The presence of true negative curvature in a state complex greatly 
increases the speed of convergence. On the other hand, a flat state complex is desirable for 
other reasons. 

It remains an important computational question to determine whether the homotopy class 
of the initial path (given, say, from a distributed algorithm) is optimal in the sense that 
its geodesic is the shortest among all homotopy classes. Interestingly, in contrast to the 
previous paragraph, this problem becomes quite a bit more difficult in the presence of 
negative curvature. From this point of view it is preferable to have a flat state complex. 



9. Discussion. Our principal contributions are: 

1. A mathematical definition of a metamorphic system which encompass many models 
currently studied and suggests seemingly unrelated (and often simpler) systems. 
Given the difficulty of building large metamorphic systems, simpler examples pos- 
sessing the same formal structure may be valuable. 

2. The state complex, whose naturality is manifested on the level of its topology (S 
can be homotopically simple) and its geometry (non-positive curvature is universal, 
mathematically helpful, and highly desirable for shape planning). 

There are drawbacks to this approach. Primary among them are the dual dilemmas that 
shape planning is inherently complex, and that there are many types of reconfiguration 
possible. A state complex approach is not meant for all systems. Indeed, it is possible to 
design degenerate metamorphic systems with little to no commutativity. Nevertheless, pay- 
ing attention to the geometry lurking behind our higher-dimensional versions of transition 
graphs leads to a nontrivial result on the optimality of path-shortening. 
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As an interesting side-note, we observe that in the present context, the optimal reconfigura- 
tion path is not always coincident with a geodesic on the configuration space. This principle 
holds in more generality. 

Finally, in this initial work, we have focused only on the first-order problem of shape plan- 
ning. We hope that the mathematical framework here suggested finds uses in more sophis- 
ticated task-planning problems for metamorphic and reconfigurable robots. 



Appendix A. Shape complexes. 

State complexes, despite their relative simplicity over the transition graph, nevertheless 
typically contain a large number of cells. For lattice-based systems, and in some other cases 
as well, a much smaller complex called the shape complex can be used to encode all of the 
local reconfigurations, without losing any essential information from the state complex. The 
idea is to exploit the symmetries of the domain. 

Consider the (lattice-based) pivoting hex system of Examplc l3.ll The underlying lattice has 
translational symmetry in two independent directions; so for instance, in the idealized case 
where the workspace is the entire (infinite) lattice with no obstacles, the state complex S 
shares these symmetries. If we take a quotient of S by the actions of these two translations, 
we obtain the shape complex Sh. (The precise definition follows shortly.) 

Observe that Sh is much smaller than S; indeed even if the workspace is infinite, the complex 
Sh is compact (provided the system consists of connected sets of, say, N cells). Yet, Sh 
carries substantially the same information as S. Whereas S keeps track of both the shape 
and the location of the aggregate, the quotient Sh ignores the location. 

Note that no information about the system is lost in the process of forming the quotient; <S 
can be completely reconstructed from Sh. Indeed, S is a covering space of Sh (with covering 
transformation group equal to the group of symmetries of the lattice) . In particular Sh and 
S share the same universal covering space, so Sh (just like S) is non-positively curved. 

Here is the definition of the shape complex; compare with Definition 14.31 

Definition A.l. The shape complex Sh of a (lattice-based) local reconfigurable system 
is the abstract cube complex whose fc-dimensional cubes are equivalence classes [U; ($ Qi )f =1 ] 
where 

1. (&ai)i=i is a fe-tuple of commuting actions of generators <j> ai ; 

2. U is some state for which all the actions ($ a< )*_j are admissible; and 

3. [Uq; (3>a 4 )*=i] = [Ui\ ($j9 4 )*=i] if and only if there is a translation T : L — ► C such 
that the list (3?^ o T) is a permutation of ($c 4 ), and Uq = U\ o T on the set 

W-U$a,(SUP(0 Q J) . 

The boundary of a /c-cube is obtained just as it is in S. 

EXAMPLE A. 2. The shape complex for the system of Example 13. II (with the generator of 
Fig. I3.1| left]l having a total of N — 3 occupied cells is shown in Figure |A~T1 There are nine 
square cells. Two pairs of edges are identified according to the matching arrows, yielding 
a Mobius strip; additionally, the three black vertices are identified to one, as are the three 
white vertices. The shape complex is therefore a Mobius strip with some of the boundary 
points identified. 
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Fig. A.l. The shape complex for a pivoting hex system consisting of 3 robotic cells. 



Macros. The most obvious advantage of the shape complex is that the discovery of a single 
shortest path in Sh simultaneously solves several state-changing problems in the complex 
S. This is particularly helpful in systems which require frequent transportation within the 
workspace (perhaps to perform various tasks in different areas). For instance, suppose the 
workspace includes several narrow corridors through which the robot must travel at various 
times. At each encounter with such a corridor, it must change to a shape which is skinny 
enough to navigate the passageway. But this is the same problem at each corridor, even 
though it occurs in different places in the state complex: the three hexagons in Example I A. 21 
can move from a triangular shape to a linear shape by the same sequence of moves, wherever 
the problem may arise in C. 

In the shape complex, we view this as a single shortest-path problem, whose solution we 
may lift to various locations in the domain. In this sense a shortest path P in Sh can be 
viewed as a macro for solving state-changing problems in S. 

A drawback to the shape complex is that it neglects the shape of the workspace W. In 
the presence of obstacles or at the edge of W, a path P in Sh may fail to lift to a path 
in S. Thus before lifting a path P one must check that the actions are admissible at the 
appropriate locations in W. 



REFERENCES 



[I] A. Abrams, Configuration spaces and braid groups of graphs. Ph.D. thesis, UC Berkeley, 2000. 

[2] A. Abrams and R. Christ, Finding topology in a factory: configuration spaces, Amer. Math. Monthly 
(109), 140-150, 2002. 

[3] A. Abrams and R. Christ, Cubical complexes for reconfigurable systems, in preparation. 

[4] M. Bridson and A. Haefliger, Metric Spaces of Nonpositive Curvature, Springer- Verlag, Berlin, 1999. 

[5] Z. Butler, S. Byrnes, and D. Rus, Distributed motion planning for modular robots with unit-compressible 

modules, in Proc. IROS 2001. 
[6] Z. Butler, K. Kotay, D. Rus, and K. Tomita, Cellular automata for decentralized control of self- 

reconfigurable robots, in Proc. IEEE ICRA Workshop on Modular Robots, 2001. 
[7] G. Chirikjian, Kinematics of a metamorphic robotic system, in Proc. IEEE ICRA, 1994. 
[8] G. Chirikjian and A. Pamecha, Bounds for self- reconfiguration of metamorphic robots, in Proc. IEEE 

ICRA, 1996. 

[9] D. Epstein ct al., Word Processing in Groups. Jones & Bartlett Publishers, Boston MA, 1992. 
[10] R. Christ, Configuration spaces and braid groups on graphs in robotics, AMS/IP Studies in Mathematics 
volume 24, 29-40, 2001. 

[II] R. Christ and D. Koditschek, Safe cooperative robot dynamics on graphs, SIAM J. Cont. &; Opt. 40(5), 

1556-1575, 2002. 

[12] M. Gromov, Hyperbolic groups, in Essays in Group Theory, MSRI Publ. 8, Springer- Verlag, 1987. 
[13] K. Kotay and D. Rus, The self-reconfiguring robotic molecule: design and control algorithms, in Proc. 
WAFR, 1998. 



COMPLEXES FOR METAMORPHIC ROBOTS 19 

[14] C. McGray and D. Rus, Self-reconfigurable molecule robots as 3-d metamorphic robots, in Proc. Intl. 

Conf. Intelligent Robots & Design, 2000. 
[15] S. Murata, H. Kurokawa, and S. Kokaji, Self-assembling machine, in Proc. IEEE ICRA, 1994. 
[16] S. Murata, H. Kurokawa, E. Yoshida, K. Tomita, and S. Kokaji, A 3-d self-reconfigurable structure, in 

Proc. IEEE ICRA, 1998. 

[17] A. Nguyen, L. Guibas, and M. Yim, Controlled module density helps reconfiguration planning, in Proc. 
WAFR, 2000. 

[18] G. Niblo and L. Reeves, The geometry of cube complexes and the complexity of their fundamental 

groups, Topology, 37(3), 621-633, 1998. 
[19] A. Pamecha, I. Ebert-Uphofl", and G. Chirikjian, Useful metric for modular robot motion planning, in 

IEEE Trans. Robotics & Automation, 13(4), 531-545, 1997. 
[20] J. Walter, J. Welch, and N. Amato, Distributed reconfiguration of metamorphic robot chains, in Proc. 

ACM Symp. on Distributed Computing, 2000. 
[21] J. Walter, E. Tsai, and N. Amato, Choosing good paths for fast distributed reconfiguration of hexagonal 

metamorphic robots. In Proc. IEEE ICRA, 2002. 
[22] M. Yim, A reconfigurable robot with many modes of locomotion, in Proc. Intl. Conf. Adv. Mechatronics, 

1993. 

[23] M. Yim, J. Lamping, E. Mao, and J. Chase, Rhombic dodecahedron shape for self-assembling robots, 

Xerox PARC Tech. Rept. P9710777, 1997. 
[24] Y. Zhang, M. Yim, J. Lamping, and E. Mao, Distributed control for 3-d shape metamorphosis, Aut. 

Robots. J., 41-56, 2001. 

[25] E. Yoshida, S. Murata, K. Tomita, H. Kurokawa, and S. Kokaji, Distributed formation control of a 
modular mechanical system, in Proc. Intl. Conf. Intelligent Robots & Sys., 1997. 































K 










• — 






• 





















1 
































• — 


► 




• 



